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GENERALIZED LINEAR BOLTZMANN EQUATIONS FOR PARTICLE 
TRANSPORT IN POLYCRYSTALS 

JENS MARKLOF AND ANDREAS STROMBERGSSON 


The linear Boltzmann equation describes the macroscopic transport of a gas of non-interacting 
point particles in low-density matter. It has wide-ranging applications, including neutron 
transport, radiative transfer, semiconductors and ocean wave scattering. Recent research 
shows that the equation fails in highly-correlated media, where the distribution of free path 
lengths is non-exponential. We investigate this phenomenon in the case of polycrystals whose 
typical grain size is comparable to the mean free path length. Our principal result is a new 
generalized linear Boltzmann equation that captures the long-range memory effects in this 
setting. A key feature is that the distribution of free path lengths has an exponential decay 
rate, as opposed to a power-law distribution observed in a single crystal. 

1. Introduction 

The Lorentz gas, introduced by Lorentz in the early 1900s [12j to model electron transport 
in metals, has become one of the most prominent objects in non-equilibrium statistical me¬ 
chanics. It describes a gas of non-interacting point particles in an infinite array of identical 
spherical scatterers. Lorentz showed, by adapting Boltzmann’s classical heuristics for the hard 
sphere gas, that in the limit of low scatterer density (Boltzmann-Grad limit) the evolution 
of a macroscopic particle cloud is described by the linear Boltzmann equation. A rigorous 
derivation of the linear Boltzmann equation from the microscopic dynamics has been given in 
the seminal papers by Gallavotti |9j, Spohn pMj and Boldrighini, Bunimovich and Sinai [2], 
under the assumption that the scatterer configuration is sufficiently disordered. 

For scatterer configurations with long-range correlations, such as crystals or quasicrystals, 
the linear Boltzmann equation fails and must be replaced by a more general transport equa¬ 
tion that takes into account additional memory effects. The failure of the linear Boltzmann 
equation was first pointed out by Golse mi for periodic scatterer configurations. We subse¬ 
quently provided a complete microscopic derivation of the correct transport equation in this 
setting mi- An important characteristic is here that the distribution of free path lengths has 
a power-law tail with diverging second moment HI m III IS] and the long-time limit of the 
transport problem is superdiffusive [[22]. Surveys of these and other recent advances on the 
microscopic justification of generalized Boltzmann equations can be found in HB Ham- 

Independent of these developments, Larsen has recently proposed a stationary generalized 
linear Boltzmann equation for homogeneous media with non-exponential path length distri¬ 
bution, non-elastic scattering and additional source terms. We refer the reader to Larsen and 
Vasques HaGUGI! and Frank and Goudon [8j for more details and applications. Prompted 
by a question of Larsen, the present paper aims to generalize our findings for periodic scatterer 
configurations mmsuiMBi to polycrystals. 

The approach of this study combines the methods of mm with equidistribution theorems 
from [2D], which were originally developed to analyse unions of incommensurable lattices. We 
will argue that, in the Boltzmann-Grad limit, the time evolution of a particle cloud is governed 
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by the transport equation 


( 1 . 1 ) 


[d t + v ■ \7 X - <%] f t (x , v, £, v+) 



ft(x, v 0 , 0, v) p o (v 0 , x , v , £, u+) cfo 0 


subject to the initial condition 


( 1 . 2 ) 


Jim f t (x,v,£,v+) = f 0 (x,v)p(x,v,£,v+), 

t— 


where fo(x,v) is the particle density in phase space at time t = 0 and p(x,v,£,v+) is a 
stationary solution of (1.1). The variables £ and v + represent the distance to the next collision 


and the velocity thereafter. By adapting our techniques for single crystals m , we will compute 
the collision kernel po(vo,x,v,£,v+) in terms of the corresponding kernel of each individual 
grain. The kernel yields the conditional probability measure 


(1.3) 


po(v 0 ,x,v,£,v + )d£dv + 


for the distribution of (£, v + ) E M>o x Sj* -1 conditional on Vo, x, v. If the grain diameters are 
sufficiently small on the scale of the mean free path length, the collision kernel has an explicit 
representation in terms of elementary functions. This yields particularly simple formulas in 
dimension d = 2. The necessity of extending the phase space has already been observed 
in the case of a single crystal Hum, finite unions [20j and in quasicrystals [T5J [21] , where 
the collision kernel is independent of ax If the scatterer configuration is disordered and no 
long-range correlations are present, the dynamics reduces in the Boltzmann-Grad limit to the 
classical linear Boltzmann equation E2H2H2H- 


The generalized linear Boltzmann equation (1.1) can be understood as the Fokker-Planck- 


Kolmorgorov equation (backward Kolomogorov equation) of the following Markovian random 
flight process: Consider a test particle travelling with constant speed along the random tra¬ 
jectory 

(1.4) x(t) = x Vt + (t — T Vt )v Vt , x(0) = xq, v(t) = v Vt , v(0) = v o , 
where 

n n 

(1.5) x n = Xq + q n , q n = 'y ^ , T n := ^ > Tj := 0, 

i = 1 i =1 


are the location, displacement and time of the nth collision, v n the velocity after the nth 
collision, and 


(1.6) vt := maxjn E Z>o : T n < t} 

is the number of collisions within time t. The above process is determined by the sequence of 
random variables (Cj, v j)jeN and (xo,vo), where (aio, woi Ci, I’l) is distributed according to 

(1.7) fo(xo, «o) p(xo, vo, fi, vi) dx 0 dv 0 d£i dv i, 

with /o now being an arbitrary probability density, and (£ n ,v n ) is distributed according to 
(1-8) po (v n —2 , X n —i, V n —i, £ n , Vn) d^ n dv n , 

conditional on {^j,Vj)^Z{ and (®o,r>o)- 

Acknowledgements. JM would like to thank Martin Frank, Kai Krycki and Edward Larsen 
for the stimulating discussions during his visit to RWTH Aachen in June 2014, and in partic¬ 
ular Edward Larsen for suggesting the problem of transport in polycrystals. We thank Dave 
Rowenhorst for providing us with the image in Fig. [lj 
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Figure 1. Grains in a sample of the /3-titanium alloy Ti-21S. The image is 
reproduced from ref. [23j| by Rowenhorst, Lewis and Spanos. 

2. The setting 

Let {Gi}i be a countable collection of non-overlapping convex open domains in M rf , and 
{Ci\i a collection of affine lattices of covolume one. We can write each such affine lattice as 
Ci = (Z d + with row vector cu,; £ and matrix Mi £ SL(d, M). We define a polylattice 

V e as the point set 

(2.1) ?> e = U {QiCeCi), 

l 

where e > 0 is a scaling parameter. We refer to Gi as a grain of V £ . An example for {Gi}i 
is a collection of convex polyhedra that tesselate i.e., U iGi = In general we will, 
however, allow gaps between grains. The standing assumption in this paper is that the number 
of grains intersecting any bounded subset of is finite. The assumption that grains are 
convex will allow us to ignore correlations of trajectories that re-enter the same grain without 
intermediate scattering. It is not difficult to extend the present analysis to include these 
effects. The assumption that all lattices have the same covolume is made solely to simplify 
the presentation and can easily be removed. Figure [IJ reproduced from [23], shows the grains 
of an actual polycrystal sample, the /3-titanium alloy TL-21S. The /3-form of titanium has a 
body-centered cubic lattice, which can be represented as the linear deformation Z 3 M of the 
cubic lattice Z 3 , where 

/ 2 1 / 3 0 0 \ 

(2.2) M = 0 2 1 / 3 0 

\ 2-2/3 2 -2 / 3 2“ 2 / 3 / 

To model the microscopic dynamics in a polycrystal, we place at each point in V e a spherical 
scatterer of radius r > 0, and consider a point particle that moves freely until it hits a sphere, 
where it is scattered, e.g. by elastic reflection (as in the classic setting of the Lorentz gas) 
or by the force of a spherically symmetric potential. We denote the position and velocity at 
time t by x(t) and v(t). Since (i) the particle speed outside the scatterers is a constant of 
motion and (ii) the scattering is elastic, we may assume without loss of generality ||u(t)|| = 1. 
The dynamics thus takes place in the unit tangent bundle T 1 (/C £!? .) where K £)r C is the 
complement of the set £>)? + V e . Here Bf denotes the open ball of radius r, centered at the 
origin. We parametrize T 1 (/C £ir ) by (x, v) £ K e ^ r x S^ -1 , where we use the convention that for 
x £ dK, t)r the vector v points away from the scatterer (so that v describes the velocity after 
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the collision). The Liouville measure on T ^(K, er ) is dv(x,v ) = dxdv, where dx = dv ol R d(®) 
and dv = dv olgd-i(w) refer to the Lebesgue measures on and S^ -1 , respectively. 

3. Free path length 

The first collision time with respect to the initial condition (x,v) £ T 1 (/C £iT .) is 

(3.1) ri(®, v) = inf{t > 0 : x + tv /C £j7 .}. 

Since all particles are moving with unit speed, we may also refer to t\(x,v) as the free path 
length. The mean free path length, i.e. the average time between collisions, is for r —> 0 
asymptotic to a - 1 e d r ~( d ^ 1 ) where if = vol/3 ^ -1 (the total scattering cross section in units of 
?’); this calculation only takes into account the time travelled inside the grains. The scaling 
limit we are interested in is when the typical grain size is of the order of the mean free path 
length. We choose (without loss of generality) e = r^ d ~ l ^ d and fix this relation for the rest of 
this paper. The mean free path length in these units is thus a -1 . 

Given (x,v) £ M. d x S]* -1 , we call the sequence (i u )i/eN the itinary of (x,v) if i v = i„(x,v) 
is the index of the iTh grain Qj v traversed by the trajectory {( x + tv,v) : t > 0}. We denote 
by £~ = £~(x,v),£+ = £f(x,v) £ [0,oo] the entry resp. exit time for Qj v . If x £ Qi x , or if 
x £ dQi x and v points towards the grain, we set £f = 0. We furthermore define the sejour 
time for each grain by £ u := £ff — £~. Note that, if a; £ Qi x and s is suffciently small so that 
x + sv £ Qi x , then 

(3.2) £f(x + sv 1 v) = £f(x 1 v) = 0, £f(x + sv,v) = £f(x,v) — s 
and, for all v > 2, 

(3.3) £^(x + sv,v) = £^(x,v) — s. 

We now consider initial data of the form (x er ,v) = (x + eq + r/3(v),v), where v £ S ^ _1 
is random, x, q £ are fixed and (3 : is some fixed continuous function. For 

x t := x + eq ^ V e (the particle is not r-close to a scatterer), the free path length Ti(x e<r , v) is 
evidently well defined for r sufficiently small. If x e £ V e (the particle is ?’-close to a scatterer), 
we assume in the following that (3 is chosen so that the ray f3(v) +M>o^ lies completely outside 
the ball B d for all v £ Sf ' 1 (thus r/3(v) + M>of lies outside B d for all r > 0). 

Let S be the commensurator of SL(d,Z) in SL(d,M). We have 

5 = {(detT)- 1/d T : T £ GL (d,Q), det T > 0}, 

cf. O Thm. 2], as well as [25,1, Sec. 7.3]. We say that matrices Mi, M 2 ,... £ SL(d, M) are 
pairwise incommensurable if MjM~ 1 ^ S for all i 7 ^ j. The pairwise incommensurability of 
Mi, M 2 ,... is equivalent to the fact that the lattices = (Z rf + u>j)Mj (i = 1,2,...) are 
pairwise incommensurable, in the sense that for any i 7 ^ j, c > 0 and u> £ the intersection 
Ci n (cCj + u>) is contained in some affine linear subspace of dimension strictly less than d. 
A natural example in the present setting would be a sequence of matrices Mi = MK{ with 
M = 1 (or M as in ( |2.2[ )) and incommensurable rotation matrices Ki £ SO(d), corresponding 
to (body-centered) cubic crystal grains with pairwise incommensurable orientation. 

The following two theorems comprise our main results for the distribution of free path 
length. The first theorem deals with generic initial data, the second when the initial position 
is near or on a scatterer, but still generic with respect to lattices in other grains. 

We will in the following use the notation 

roc r£ 

(3.4) D&(£) = / ${ri)drj = 1 — / ^(ri)drj 

Jo 

for the complementary distribution function of the probability density <L. 

In the following we consider lattices Ci = e~ 1 x + (Z rf + ujffMi with an additional shift by 
e~ l x. This looks artificial but is necessary for all subsequent statements to hold. [The problem 
becomes easier if we assume that u>i are independent random variables uniformly distributed 
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in the torus M d /Z rf . In this case it is fine to use Li = (Z d + Ui)Mi. All of the statements below 
will also hold in this case.] 

Theorem 1. Fix x G and, for all i G N, let Li = e~ l x + (Z d + Ui)Mi with u>i G M. d , 
and Mi G SL(d, M) pairwise incommensurable. Fix q G so that Ui — qM- _1 ^ /or all i. 
If (v)i/eN the itinary of ( x , v), then, for any Borel probability measure X on S ^ -1 and any 


e > 0 , 



(3.5) 

POO /• 

lim \({v G S^ 1 : Ti(® e , r ,v) >/})=/ / 

Jt Jsf- 

\l/( x ,u, q) d\(v) dq 

with 



(3-6) 

w(x ,„, e) = ((n?R.M*K-a 

if If <(< t-i 

otherwise, 


where <!>(£) is the limit probability density of the free path length in the case of a single lattice 
and for generic inital data, see [HI Eq. (1.21)]. 


By [HI Eq. (1.23)] we have, 

(3-7) ^) = a-^ + 0(e), 


where ((d) is the Riemann zeta function and the remainder is non-negative. In dimension 
d = 2 the error term in fact vanishes identically for £ sufficiently small; indeed, for 0 < £ < ^ 
we have [H Theorem 2] 

24 

(3-8) $ (e) = 2 --£ 

7 T z 

and hence 

(3.9) D*(£) = l-2£ + ^/ 2 . 

In dimension d = 3 we have [TH] Corollary (1.6)] for 0 < ( < | 


(3.10) 


$(£) = 7T - 


7T 


C (3) 


£ + 


37T 2 + 16 c 2 
2vrC(3) ^ 


and so 
(3.11) 




1 - ir( + 


7r 


2C(3) 


37T 2 + 16 3 
6ttC(3) 


This means that the limit distribution T(a;, v, () is completely explicit, if the diameter of each 
grain is bounded above by 2 in dimension d = 2 resp. | in dimension d = 3. 

Let us now turn to initial data near a scatterer. Let us fix a map K : S^ 1 —> SO (d) such 
that vK(v) = e.\ for all v G where e\ := (1,0,..., 0). We assume that K is smooth when 
restricted to S^ 1 minus one point (see m footnote 3, p. 1968] for an explicit construction). 
We denote by x± the orthogonal projection of x G M. d onto the hyperplane perpendicular to 
ei- 


Theorem 2. Fix x G Qj for some j G N, and, for all i G N, let Li = e -1 a? + (Z d + uji)Mi with 
u>, G and Mi G SL(d, M) pairwise incommensurable. Fix q G such that x e = x + eq G 
eLj and such that oji — qMf 1 ^ Q d for all i 7 ^ j. If (i v )u &N is the itinary of (x, v), then, for 
any Borel probability measure A on S ^ 1 and any ( > 0, 


(3.12) lirn A({u G Sf 1 : Ti(x e<r , v) > £}) 



f ^>o(x,v,q,(P(v)K(v))±)dX(v)dq 

J Sf - 1 
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with ^o(*, v, £, w) defined for any x G U jQj by 

'$o(£,w) 

(3.13) T 0 (x,v,£,w)=< 


ifO<£<£+ 


^i^)(np2^(4))^«-a ^ c<e<e^> 2 ; 

0 otherwise, 

where $o(£) , rO is the corresponding limit probability density in the case of a single lattice, and 
$a(r},w)dri. 

The single-lattice density is given by 4>o(f;, w ) = f B d-i < I > o(^, w, z ) dz with <I>o(£, w, z) as in 

[TH1 Sect. 1.1] (cf. also Thm. [4] below). We extend the definition of To(®, v, £, w) to all x G M d 
as follows. Given a grain Qj and (x,v) with x G dQj, we say v is pointing inwards if there 
exists some eo > 0 such that {x + ev : 0 < e < eo} C Qj. Let 

(3.14) Hj := {(x,v) G dQj x Sf _1 : v is pointing inwards} 
and 

(3.15) g 3 := (Gj x S d l - y )UU J . 

We now extend the definition of To(a;,u,^,m) to all x G M d , £ > 0, by setting 

lim e _> 0 + 4'o(* + ew, v,f — e, w) if (x, v ) G U jHj 
0 if (x,v) £l)jGj. 


(3.16) d'o (x,v,£,w) = 

Let us furthermore define 

(3.17) t(x,v) = 


1 if (x, v) G UiQi, 

0 otherwise, 

and the differential operator V by (assume £ > 0) 

(3.18) V^f(x, v, £) = lim e~ 1 [T(a; + ev,v,£ — e) — \H(x, u, £)]. 

e—>0+ 

Note that we have T>^{x,v,f) = [u • V* — e?f] ^(aq v, £) wherever the right-hand side is well 
defined (which is the case on a set of full measure). In the case of a single lattice, we have 
Eq. (1.21)] 


(3.19) 


r 

m r l 




4>o {Vi w ) dw dr\. 


In the case of a polylattice, (3.19) generalizes to 

{V^{x,v,i) = f B d-i y 0 (x,v,£,w)dw (f>0) 

[^(a:, u, 0) = al(x,v). 

This relation follows from (3.6), (3.13) and ( |3. 19 ) in view of the relations (3.2), (3.3). 

4. The transition kernel 


To go beyond the distribution of free path length, and towards a full understanding of the 
particle dynamics in the Boltzmann-Grad limit, we need to refine the results of the previous 
section and consider the joint distribution of the free path length and the precise location on 
the scatterer where the particle hits. 

Given initial data ( x,v ), we denote the position of impact on the first scatterer by 
( 4 . 1 ) X\(x, v) := x + T\(x, v)v. 

Given the scatterer location y G V e , we have x\(x,v) G S +y and therefore there is 
a unique point w\(x,v) G S^ _1 such that x\{x,v) = rw\{x,v) + y. It is evident that 
—Wi(x,v)K(v) G S', d 1 , with the hemisphere 1 = {v = (ui,..., vfi) G S^ -1 : v\ > 0}. 
The impact parameter of the first collision is h = (w\(x, v)K(v))±. 
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As in Section^ we will use the initial data (® e>r , v) = (x + eq + r/3(v),v), where i)£Sf 1 
is random, x, q E are fixed and (3 : is some fixed continuous function. 

We again have two theorems, the first for generic initial data, the second when the initial 
position is near or on a scatterer, but still generic with respect to lattices in other grains. 
Theorems [l] resp. [^follow from Theorems [i] resp. [4] below by taking the test set if = S ) d 1 . 

Theorem 3. Fix x E and, for all i 6 N ; let Li = e _1 x + (Z rf + ujj)Mj with Ui E 
and Mi E SL(d,M) pairwise incommensurable. Fix q E so that u — qMff 1 (f <Q) d /or 
all i. If (iv) U £T$ is the itinary of (x,v), then for any Borel probability measure A on Sf -1 
absolutely continuous with respect to vol g d-i, any subset if C S} 1 with vol g d-i(3il) = 0, and 
any 0 < a < b, we have 


(4.2) lim A({u E 1 : T\ E [a, b), —WiK(v) Ell}) 



where 


(4.3) 


^(x,v,C,w) 


(n if if <f<i+ 

0 otherwise. 


Theorem 4. Fix x E Qj for some j E N, and, for all i E N, let Li = e -1 ® + (Z d + u,;)Mj with 
cji E and Mj E SL(d, M) pairwise incommensurable. Fix q E such that x e = x + eg E 
e£j and sued that cji — qMff 1 Q d for all i / j. If (i v )ve N ide itinary of (x,v), then 
for any Borel probability measure A on absolutely continuous with respect to vol g d-i, any 

subset it C 1 with vol g d-i(5il) = 0, and any 0 < a < b, we have 


(4.4) lim A({u E Sf 1 : T\ E [a, b), —W\K{v) E it}) 

r—»0 v 1 ’' 



a Jil ± J S“ 


To ( x , v, £, w, (f3(v)K(v))±) d\(v) dw df 


with 


(4.5) 3> 0 {x,v,£,,w,z) = < 


$o(£,w,z) if 0<f<£+ 

Wi,z)( */c<e<# (»>v, 

0 otherwise, 

where 4>o(/i w, z) is the transition kernel for a single lattice, cf. [181 Sect. 1.1]. 

As above, we extend the definition of To(®, v, £, w, z) to all x E M. d by setting 


(4.6) T 0 (®,v,£, w,z) = 


lim e _ s . 0+ T 0 (® + ev, v, f — e, w, z) if (®, v ) E U jUj 
0 if (x,v) £ LijQj. 


We refer the reader to [16l ITBlj for a detailed study of 4>o(£, w, z), T(£, iu) and <f>o(£j t^), 
which are related via [T7| Eq. (6.67)], 


(4.7) 

and 

(4.8) 


$(£,w) = 




4>o(r/, w, z) dz dq 


®o(£,w) = / T 0 (f,w,z)dz. 

JB f " 1 
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We have in particular 
(4.9) 


Eq. (1.18)], 


1 - 2 
C (d) 


< $o(£,w,z) < 


((d) 


that is, •ho ((,w,z) = ((d) 1 + 0(£), and jTS] Eq. (1.19)] 


(4.10) 


a 




where the remainder term is everywhere non-negative, and the implied constant is independent 
of w. As for the free path lengths, we have explicit expressions for these transition kernels in 
dimensions two and three, which will be discussed in Sections [7] and [8] 

The generalization of (4.7) is 

(V^(x,v,( } w) = f B d -1 T 0 (x,v,(,w,z)dz 

|T(a;, u, 0, w) = t(x,v). 

Its proof is analogous to ( |3.20| ) . 

The single-crystal transition kernel satisfies the following invariance properties PI: 

(4.12) 4>o(£, z, w) = 4>o(£, w , z ), 
and for all R 6 O (d — 1), 

(4.13) $o(£, wR, zR) = $o(£, w, z ), 

(4.14) $(£, wR) = 4>(£, w), 4»o(?, wR) = <f> 0 (£, to). 

These relations imply 

(4.15) ^ 0 (x + (v,-v,(,z,w) = ifr 0 (x,v,(,w,z), 
and for all R € 0(d — 1), 

(4.16) ^o(®, v, (, wR , zR) = ^ 0 (®, v, (, w, z), 

(4.17) &(x,v,(,wR) = &(x,v,(,w), ^ 0 (x,v,(,wR) = ^f Q (x,v,(,w). 

5. Random point processes and the proof of Theorems SHI 


We follow the same strategy as in [16] but use the refined equidistribution theorems for 
several lattices from [20]. Recall (2.1), namely V £ = Uj (Q t n e£j), with affine lattices £,, = 


e l x 


+ (Z d + uJi)Mi. Let 


(5.1) 


A f = 


IQ e -l/(d-l)l d _ i 


E SL(d,R). 


The idea is to consider the sequence of random point processes (with e = A d L/ rf ) 

(5.2) & £ (x £ir ,v) := e~ 1 ([V t - (x + eg)] \ {0} - r(3(v))K(v)A e 

(where v is distributed according to A) and prove convergence, in finite-dimensional distribu¬ 
tion, to a random point process as e —> 0. Note that the removal of the origin in (5.2) has an 
effect only when Q; := ui — qMf 1 E Z rf for some i\ in fact we have 

(5.3) <d £ (x £jr , v) = |J - cc £ , r ) n ((Z d + ati \ {0 })M t - e 1/(d ~ 1] (3(v))\l<(v)A e . 

i ' ' 

Set G = ASL(d, M), T = ASL(d, Z), and let g be the unique G-invariant probability measure 
on T\G. We let O be the infinite product space R = fl* T\G (one factor for each Gi) and let 
uj be the corresponding product measure U; l 1 - Let us define, for any (gi) G f l and v E S^ _1 , 

(5.4) &(x,v,( 9i )) :=U[(((g i - a i)A(u) n Me 1 ) x R^ 1 ) nZ^' 











PARTICLE TRANSPORT IN POLYCRYSTALS 


9 


Theorems [3] and [4] (and thus Theorems [T] and [2]) follow from the next two theorems by the 
same steps as in [20, Sections 6 and 9]. 

Theorem 5. Fix x £ and, for all i G N, let Li = e~ l x + (Z rf + cjj)Mj with Ui £ 
and Mi £ SL(d, M) pairwise incommensurable. Fix q £ M. d so that uji — qMff 1 (f Q d . Then 
for any Borel probability measure A on S^ 1 absolutely continuous with respect to volgd-i, any 
bounded sets £>i,... ,Bf. C with boundary of measure zero, and m i,..., £ Z>o, 


(5.5) lim A({t> 6 Sf 1 : #(O e (x^ r ,v) n Bi) = mi (Vi = 1,, k)}) 


[ w({(si)En : #(G(x,v,(gi))r\Bi) 

Jsf - 1 


mi (Vi = 1, • • .,&)}) d\(v). 


Now set Go = SL(d,M), To = SL(d, Z), and let no be the unique Go-invariant probability 
measure on To\Go; then let ir-F = (To\Go) x fli/j F\G and let w be the corresponding 
product measure no X Y\,^j l 1 - Finally let us define, for any (gf) £ and v £ S^ _1 , 


eW(x,v, ( 9 i)) := [(((£,■ - x)K(v) nRd) x M d_1 ) D {(Z d \{0}) 9j - (f3(v)K(v)) ± ) 


(5.6) 


u U[(((^ -*)-K'(«)nRei) xR^nZ^ji 


Theorem 6. Fix x £ Qj for some j £ N, and, for all i £ N, let Li = e~ x x + (fL d + ojfjMi with 
u> t £ and Mi £ SL(d, M) pairwise incommensurable. Fix q £ such that x e = x + eq £ 
eLj and such that u>i — qMff 1 (f Q d for all i ^ j. Then for any Borel probability measure 
A on absolutely continuous with respect to vol g d-i, any bounded B\,...,Bk C with 
boundary of measure zero, and m \,..., mk £ Z>o, 


(5.7) limA({weSf l : #(0 e (cc e , r , v) n Bi) = mi (Vi = 1, ..., fc)}) 

= / w({( 5 i) <E : #{&( J \x,v,{gi))nBi) =m t (Vi = 1,... , k)}) dX(v). 

J sf _1 

Theorems [5] and [6] are implied by [201 Theorem 10] by the same arguments as in jl6l Section 

6 ], 


6. Tail estimates 


We will now show that, unlike the case of single crystals, the distribution of free path 
lengths, as well as the transition kernels, decay exponentially for large f. This observation 
relies on the following bound. 


Lemma 7. For £ > 0, 

(6.1) _D$(£) < max (e“^,e a - ). 

Proof. Since 1 — x < e~ x for 0 < x < 1, we have 

(6.2) D$(£) < e ~fo^ dr > 


< e 


~ fo max ( CT - dr l 


where the second inequality follows from the positivity of the error term in (3.7). 


□ 


By grain diameter we mean in the following the largest distance between any two points in 
a single grain. We define the gap function, gap(a;, v, £), to be the total length of the trajectory 
{x + tv : 0 < t < £} that is outside U jQj. Note that g&p(x,v,f) < £, and furthermore 

(6.3) gap(® + sv, v, f) = gap(®, v,£ + s)- gap(*, v, s ) 

for all s > 0. 
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Figure 2. The functions F(t) and G(w). 


Proposition 8. Assume that all grain diameters are uniformly bounded. Then there are 
constants C, 7 > 0 such that for all x, v, £, w, z 

(6.4) ^ 0 (®, v , £, w, z ) < Ce- 7( «- gap(a: ’ ,, ’ ?)) . 

The same bound holds for '!'(*, v, £, w), 4'o(ai, v, £, w) and T(ai, v, £)• 

Proof. If the grain diameters are bounded above by £ > 0, we have li < £ for all i, and hence 
in view of Lemma [7j 

(6.5) D*(£i) < e“ 7 * 


for all i, where 7 = min(^. The desired bound now follows from (4.5). 


□ 


Therefore, if gap(:c,t;,£) < 5 £ for some 5 6 [0,1), we have exponential decay in (6.4) with 
rate 7(1 — 5). 


7. Explicit formulas for the transition kernel in dimension d = 2 
In dimension d = 2 we have the following explicit formula for the transition probability [15] : 


(7.1) 
with 

(7.2) 


, lt , 6 / i 1 - max(|u>|, \z\) - 1 

$0 (£,w,z) = 7T 1+- 


7T 


T(x) = < 


| w + z\ 

0 if x < 0 
x if 0 < x < 1 
1 if 1 < x, 


The same formula was also found independently by Caglioti and Golse [7j and by Bykovskii 
and Ustinov [5], using different methods based on continued fractions. In particular, for all 


(7.3) 


®o(£,w,z) = 


7T 


which is thus independent of w,z. We have furthermore pa, again for all £ < |, 


(7.4) 


12 

$0 (£,w) = - 5 , 


12 


7T 


^(e, t«) = 1 — o e- 


7T 


Recall that in dimension d = 2, the value \ is precisely the mean free path length. 
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8. Explicit formulas for the transition kernel in dimension d = 3 
The results in this section are proved in |18j. For 0 < t < 1 set 
(8.1) F(t) = 7r — arccos(i) + t\J 1 — t 2 = Area({(xi, x-z) E B\ : xi<t}), 

see Fig. [ 2 ] Then, for 0 < £ < 


( 8 . 2 ) 

and 

(8.3) 

where G : [0,1] 


$ 0 (f,to,«) = C(3) 


-1 




TT 


$(£,*>) = 1 - + 


6 


C(3) 7r2 C(3) 


G(IMI)£ 2 , 


M>o is the function 

rl—w 


r dr, 


ft-w ri+w 2 1 

(8.4) G(w) = 7T / F(\r)rdr+ / F(fr) arccosl---) 

Jo " Ji-w ' 2 ivr ) 

cf. Fig. J 2 J The function G(w) is continuous and strictly increasing, and satishes (7(0) = 
^(4tt+3v / 3) and = ^2 + L 

9. The transport equation 

In order to prove that the dynamics of a test particle converges, in the Boltzmann-Grad 
limit r -A 0, to a random flight process (x(t),v(t)) defined in (1.4)—(1.6), we require technical 
refinements of Theorems [3] and [4j where the convergence is uniform over a certain class of 
A. This argument follows the strategy developed in for single crystals. We will here not 
attempt to prove these uniform versions, but move straight to the description of the limit 
process which is determined by the transition kernel of Theorem [4| 

As in the case of a single crystal EE!, the limiting random flight process becomes Markovian 
on an extended phase space, where the additional variables are 

( 9 . 1 ) m = T Vt+ \ — t E M>o (distance to the next collision) 
and 

(9.2) v +(t) = V v t +1 £ S^ -1 (velocity after the next collision). 

The continuous time Markov process H(i) = (x(t),v(t),£(t),v+(t)) is determined by the initial 
distribution fo(x, v, £, v + ) and the collision kernel po(vj-i,Xj, Vj,£j + i,Vj + i) which yields the 
probability that the (j + l)st collision is at distance from the jth collision, with subsequent 

velocity Vj + 1 , given that the jth collision takes place at Xj and the particle’s velocities before 
and after this collision are Vj-i and Vj. The collision kernel po(vq, x, v, £, n+) is related to 
the transition kernel To (x,v,£,w, z) of the previous sections by 

(9.3) p 0 (v 0 ,x,v,£,v+) = '& 0 (x,v,£,b,-s)<T(v,v + ) 

where a(v,v+) is the differential cross section, s = s(v,v 0 ), b = b(v,v+) are the exit and 
impact parameters of the previous resp. next scattering event. 

The density ft(x, v, £, v + ) of the process at time t > 0 is given by 


(9.4) 


/ ft(x, v, £, u + ) dx dv d£dv + = P(F(t) e A) 

Ja 


for suitable test sets A. Let us write 
(9.5) f t (x,v,£,v+) = 

where f^ n \x,v,^,v + ) is the density of particles that have collided precisely n times in the 
time interval [0, t] . Then 


OO 

J2fl n \x,v,t,v + ) 

n=0 


( 9 . 6 ) 


fi°\x,v,^,v + ) = fo(x - tv,v,£ + t,v+), 
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and for n > 1 


(9.7) fi n) {x,v,^,v + ) = f fo{x 0 ,v 0 ^i,Vi) 

JT n <t 


X n Xj, Vj, £j + i, Vj+i) din dv n -1 ■ • • di\ dv 0 , 

j=i 

with w n+ i = v + , v n = v, i n+1 = i + t - T n , x 0 = x - q n - (t - T n )v, and with x j} q n , T n 
as in S For general densities fo(x,v,i,v+), relations ( |9.5[ )-( [9?7| ) define a family of linear 
operators (for t > 0 ) 

Kt n) fo(x,v,i,v + ) := rf n \x,v,i,v + ) 


(9.8) 
and 

(9.9) 

One can show that 

(9.10) 


K t f 0 (x,v,£,v +) := f t (x,v,i,v + ). 


E rAn-m) K (m) _ K (n) 
m =0 


which in turn implies K t] +t . 2 = K tl K t2 , i.e., the operators K t form a semigroup (reflecting the 
fact that S(i) is Markovian). The proof of this is analogous to the computation in |17l Sect. 
6.2], Hence, for h > 0 we have ft+h = Khft■ Since the probability of having more than one 
collision in a small time interval is negligible, we have for small h (cf. [171 Sect. 6.2]) 

(9.11) f t+h (x,v,£,v+) = K h f t (x,v,£,v+) = f t (x, v,£, v + )+K^ ] f t (x, v, £, v + )+0(h 2 ). 


Explicitly, we have by ( |9.6[ ), ( |9.7[ ), 

(9.12) f t+h {x, v , i, v+) = f t {x - hv, v,i + h, v + ) 


+ 


' 0 j sf- 1 


ft(x - iiv 0 - (h - ii)v, v 0 , i 1 ,v) 


x p 0 (v 0 , x-(h- ii)v, v,i + h- ii,v+) dv 0 d£i + 0(h 2 ). 

Dividing this expression by h and taking the limit h —> 0, we obtain the Fokker-Planck- 
Kolmogorov equation (or Kolmogorov backward equation) of the Markov process E(t), 


(9.13) 


where 


= / f t (x,vo,0,v)po(vo,x,v,i,v + )dv 0 , 

lad -1 


(9.14) Vf t (x,v,i,v + ) = lime 1 [/*+«(* + ev, v, £ - e, v+) - f t (x, v, f,«+)]. 

£-> 0 + 

As for T>, we observe that V = dt + v ■ \7 X — d £ at any point where the latter operator is well 
defined (which is the case for a full measure set). The physically relevant initial condition is 

(9.15) lim f t (x,v,i,v + ) = f 0 {x,v,i,v + ) = f 0 (x,v)p(x,v,£,v+) 
t—>o 

with 

(9.16) p(x,v,i,v + ) := V(x,v,£,b) <r(v,v+). 

The original phase-space density is recovered via projection, 

(9.17) 


/•OG n 

ft{x,v)= / f t (x,v,£,v+)dv+d£. 

Jo J s? _1 


/ o J s“ 

Note that ( 4.11[ ) implies that ft(x, v , £, v + ) = p(x, v , £, u + ) is the stationary solution of ( |9.13[ ), 
corresponding to fo(x,v) = 1. Uniqueness in the Cauchy problem (9.13) (9.15) follows from 
standard arguments, cf. m Section 6.3]. 
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The generalized linear Boltzmann equation (9.13) will also hold for other grainy materi¬ 
als, provided different grains are uncorrelated to guarantee the factorization of the individual 
grain-distribution functions in (4.5). We will discuss the simplest example, grains of a disor¬ 
dered medium, in Section 11 Note that it is not necessary that the transition probabilities 
<l>o(£, w i z ) i n each grain are identical—the modifications in ( |4.5| ) are straightforward: replace 
4>o(£, w, z), $(£,«>), etc. by the grain-dependent $q (£, w,z), w), etc. throughout. 

10. A SIMPLIFIED KERNEL 

Let us assume now that the transition kernel is given by a function 

(10.1) ^o(*,v,0 := v^o(x,v,£,w,z) 

that is independent of w,z. As we have seen in Section [7j this holds in the two-dimensional 
setting provided the grain size is less than the mean free path length. Then of course also 
'L(a:, v, £, w) is independent of w and related to the distribution of free path lengths by 

(10.2) T(a:, v, £) = a^f(x, v, £, w). 


Relation (4.11) becomes (3.20); that is 

( V^/(x, v, £) = a^ 0 (x,v,^) 


(10.3) 

The ansatz 

(10.4) 


1 *(*, v, 0) = a l(x, v ). 
ft(x, V, £, v+) = <j~ 1 g t (x, v, f) <j(v, v + ) 


reduces equation (9.13) to 

(10.5) Vg t {x,v,£) = (j- 1 Vo(x,v,€) / g t (x, v 0l 0) cr(v 0 , v) dv 0 

J s ^ 1 

with initial condition 

( 10 . 6 ) Kmgt(x,v,€) = f 0 (x,v)V(x,v,£). 

t-> o 

The stationary solution of ( 10. 5| ) corresponding to fo(x,v) = 1 is gt(x,v,£) = \k(:r, u, £). 

11. Disordered grains 

It is instructive to contrast the case of crystal grains discussed above with grains consisting 
of a disordered medium. We model the medium by scatterers centred at a fixed realisation of 
a Poisson point process £ Poisson with intensity 1, rescale by e and intersect with the grains as 


in (2.1) to produce the point set 

( 11 . 1 ) 


v e = (Qi n eCp oisson) • 


If there are no gaps between the grains, V e is precisely a fixed realisation of a Poisson process 
with intensity 1, for which the convergence to the linear Boltzmann equation has been estab¬ 
lished in [2], There do not seem to be any technical obstructions in extending these results 
to the setting with gaps. In particular, the above results for crystal grains remain valid, if we 
replace the relevant single-crystal distributions by their disordered counterparts (cf. [13], [T3]): 
For the transition kernels, we have 

(11.2) $o(f, w, z ) = e“ CT? , $(£, w) = e“^, 
and for the distribution of free path lengths 

(11.3) $o {£,w) = ae~ a ^ $(£) = ue _<7? , D$(£) = e~^. 


Thus in the case of a disordered granular medium, the formulas (3.6), (3.13), (4.3) and (4.5) 
become 

(11.4) ^(s, v,£,w) = <7 ^(x, v, £) = e _IT ^ _gap ( a: ’' u ’^)l(cc + £v, v), 
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and 

(11.5) l I / o(a:, v, £, w , z) = a _1 ^o(®, v, £, w>) = e _cr ^ _gap ^ 3: ’' !, ^^l(a;, v)l(* + £v, v ). 

The kernel \I/o(#jV,£) := <r d'o(®, w, £, 2) in (11.5) is evidently independent of and 


thus we are in the setting of Section [TO] We have 

(11.6) ^o(®, = ^(*,v,£)l(*,v), 

and the ansatz g t (x,v,£) = ft(x,v)^(x,v,^) reduces the generalised Boltzmann equation 
( 10. 5| ) to the classical density-dependent linear Boltzmann equation 

(11.7) [d t + v ■ V x ]f t (x, v) = t(x, v) / [ft(x,v 0 ) - ft(x,v)\ a(v 0 ,v)dv 0 

/ od-1 


with initial condition 

( 11 . 8 ) 


s“- 


]im ft(x,v) = f Q (x,v) 
t—>0 
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